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We analyze a proposed scheme of detection of phase-shift-keyed signals 
using an injection-locked oscillator the bandwidth of which is much less 
than the modulation rate. The output of the oscillator is a carrier with 
essentially all of its modulation removed. We analyze the effect of noise and 
signal modulation on the phase of this reference tone and compute its effect 
on the probability of detection error. If a suitable encoder and decoder are 
used for the transmitted signal, this technique can provide nearly ideal 
coherent demodulation. 

I. INTRODUCTION 

The two generally recognized methods of detection of phase-shift- 
keyed (PSK) signals are coherent detection and differential detection. 
Coherent detection has been shown to be optimum in the presence of 
Gaussian noise, 1 but, due to the difficulty of storing an absolute phase 
reference at the receiver, it is seldom used in practice. 

In a recent paper, 2 B. Glance showed that an injection-locked oscil- 
lator, the locking bandwidth of which is much less than the modulation 
rate can, under certain conditions, be used to derive a phase reference 
from the input signal itself. This is actually a form of a quadrature 
reference system, where the phase of one quadrature remains essentially 
unkeyed, and is used to provide the reference tone. 3,4 In this paper, we 
examine the effectiveness of this scheme for a two-phase PSK system 
where the modulation is a random digital signal and additive Gaussian 
noise is present. We derive an expression for the probability distribution 
of the reference phase, and from this calculate the average probability of 
a detection error. We find that if the modulation rate is much greater 
than the bandwidth of the oscillator, and if a suitable encoder and 
decoder are used, the method very nearly approaches the ideal perform- 
ance of coherent detection. 
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II. LOCKING EQUATION ANALYSIS — ZERO ORDER SOLUTION 

A portion of the received signal is used as the input to the injection- 
locked oscillator. This signal may be represented as 

x(t) = V2 A cos M + 0(!)\ + »(0i C 1 ) 

where A is the signal power, co is the carrier frequency, and 0(0 is the 
phase modulation. The received signal is assumed to be contaminated by 
additive white Gaussian noise, n{t), with double-sided spectral density 
N /2. If co is sufficiently close to the natural oscillator frequency, w , 
locking will occur and the output of the oscillator will be 

y(t) = V2 B cos M + 0(0 - 0(0], (2) 

where B can be assumed to be constant. 5 <j>(t) is the phase difference 
between the input and output signals of the oscillator. In the case of 
interest, the total phase modulation of the oscillator output, rj(f) = 
g(t) _ ^(t), is small, and y(t) is used as the phase reference in the coherent 
detection of the remaining portion of the received signal. 

In the absence of noise, the phases of the input and output signals of 
the oscillator are related by the well-known locking equation 6 

^f + A sin 4>(t) = co - co + —gr » W 

where 2A is the locking bandwidth of the oscillator. This equation takes 
the same form as that for a first order phase-locked loop. 7 The effect of 
the noise at the input has been analyzed by Viterbi. 8 The effect is to add 
an additional term to eq. (3), 

where n'(t) has the same statistics as n(t). We rewrite eq. (4) in terms of 
v(t). 

^ + A sin [„(0 - 6(t)] = coo - o> + J n'(t). (5) 

For the case of zero input phase modulation, i.e., 6(t) = 0, eq. (5) 
becomes 

^ + A Bin „(<) - Wo - « + J n'(t). (6) 

Using Fokker-Plank techniques, Viterbi derived from this equation 
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p(v), the steady-state probability density of 77. For to = o> the solution is 

pW - S)' (7) 

where 7 (- ) is the zeroth order modified Bessel function of the first kind> 
and a 2 = 4:A 2 /N A is the signal-to-noise ratio in the bandwidth of the 
oscillator. For a » 1, this distribution for 77 small is nearly Gaussian 
with mean zero and standard deviation l/a. For u ^ w the distribution 
becomes centered about the point /3 = sin -1 (« — <*>)/A. In the case, 
a » 1 and |(w () — co)/A | « 1, ^ « (w - W )/A and for 77 small, the 
distribution is very nearly equal to 

We now consider the case 0(0 9^ 0. In a binary PSK signal, 0(0 is a 
waveform of the form 

0(0= Y,a nV {t-nT), (9) 

where a„ = ±1, and p(t) is assumed to be a pulse which is nonzero only 
over the range < t < T. For the moment we assume zero noise. The 
resulting equation is 

^p + A sin [7,(0 - 0(0] =0,0-0; (10) 

which we rewrite as 



dy(t) 
dt 



+ A sin t;(0 cos 0(0 — A cos tj(0 sin 0(0 = o> — o>. (11) 



Our technique for the solution of this case will suggest a method of 
handling the stochastic problem when the noise term is reintroduced. 

The exact solution to eq. (11) is difficult or impossible to obtain for 
general 0(0- Let us therefore take advantage of the fact that T <3C 1/A to 
derive a differential equation, the solution of which approximates that 
of eq. (11). 

We assume that there exists an interval of length r with the property 
that T <<C t <3C 1/A, and we take the average of eq. (11) over t. Letting ( ) 
denote this averaging operation, i.e., 

<n(0) -- f" ' v(u)du (12) 

T J l-(r/2) 
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there results 

^P + A(sin ij(0 cos d(t)) - A(cos v (t) sin d(t)) - <*„ - «. (13) 

Choosing t « 1/A insures that r/(0 is very nearly constant over the 
interval of averaging, for, from eq. (10), we have \drj{t)/dt\ g A + 
|w - <a\ g 2A. Consequently, ij(<) » <i?(0) and the quantities sin q(t) 
and cos 77(0 may be taken outside the averaging operator. 

^4^ + A sin (v(t))(coa 0(t)) - A cos (v(t))(am 0(f)) « w - w. (14) 

On the other hand, the choice t :» T insures that there will be many data 
pulses over the interval of averaging. Thus the quantity (cos 0(0) is very 
nearly constant and is equal to 



C = (cos 0(0) = I f cos p(0 d*. 



(15) 



(Table I lists the value of this quantity for three important pulse shapes.) 
Since we usually have no control over the transmitted message, the 
quantity (sin 0(0) will not, in general, be time-independent. However, a 
scheme has been suggested by C. L. Ruthroff and W. F. Bodtmann 9 in 
which, through the use of a simple encoder and decoder, this quantity 
can be made very nearly equal to zero. 

Table I— Values of C and 0- for Three Pulse Shapes 
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The encoder is a device which stores a block of N bits of the message, 
and counts the net difference of +l's and — l's contained in the block. 
It also keeps a separate count of the net difference in +l's and —l's 
which have been sent over the entire past history of the message. The 
entire block is transmitted either with normal polarity, or with reversed 
polarity, in such a way as to cause the accumulated count to come as 
close as possible to zero. Preceding each block is a single "code" bit which 
specifies the polarity of that block. (The code bits are included in the 
counts.) The decoder decodes the message in an obvious manner. 

We assume that 0(0 is the output signal of such an encoder. Assuming 
t 2> NT, we have 

(sin 0(0) « 0, (16) 

so that our approximating differential equation becomes 



dy (t) 
dt 



-f- AC sin rj o (0 = co — cj. 



(17) 



We call the solution to this equation the "zero-order approximation to 
ri(t)". 

We note that this equation has the identical form as eq. (10) for the 
case of zero modulation. Thus the "zero-order effect" of the modulation 
is to reduce the effective value of the locking bandwidth by a factor C. 

In Figs. 1 and 2, we demonstrate the above results. In both figures, 
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Fig. 1 — »j(0 with square-wave modulation. 
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Fig. 2 — 77(0 with random modulation. 

A = tt/807 1 . The lower curve of Fig. 1 shows the behavior of rj(t) versus 
time for the case of zero modulation, i.e., 0(0 = 0. The saw-tooth wave 
shows the behavior of tj(0> that is, the exact solution to eq. (10), for the 
case where 0(0 is a square wave with amplitude 3ir/8 and period 4T. 
The smooth curve drawn over the saw-tooth wave shows the zero-order 
approximation obtained from eq. (17). As can be seen, the main effect of 
the modulation has been an increase in the decay time of the resulting 
curve. This is a result of the decrease in the effective locking bandwidth 
caused by the modulation, as predicted above. 

In addition to this, the true curve has a "wiggle" which seems to 
increase in size as the curve approaches zero, and which reaches a 
maximum magnitude of about 0.02 X r/2. 

Figure 2 shows the true curve and the zero-order approximation where 
the phase modulation is a random binary signal with Prob (+1) = 
Prob (—1) = \ which has then been passed through an encoder with 
N = 5; and a rectangular pulse shape of amplitude 3tt/8 is used. We 
note the curve follows the same overall path as in Fig. 1. The "wiggles" 
are now random in character; we note that their amplitude again appears 
to grow as the curve approaches zero. 

III. FIRST ORDER SOLUTION 

To better understand the behavior of these "wiggles" we now derive a 
correction term, 1^(0. which, when added to y (Q, improves the accuracy 
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of our solution. Wo call 77,, (0 + 77 ,(0 the "first-order approximation to 

We define j,'(t) = v(0 - v,.(t). Subtracting eq. (17) from eq. (10), 
there results 

^- + A sin r, f (t) cos [, o (0 - d(t)] 

= AC sin ijo(0 - A cos rj'(t) sin [ij o (0 - 0(0] • (18) 

If the assumptions we have made above are valid, the error in the zero- 
order approximation will be small, |tj'(0| « 1, and we may linearize eq. 

(18). 

^ + A„'(0 cos [, o (0 - 0(0] 

= AC sin r, (t) - A sin [ Vo (t) - 6(1)]. (19) 

Since IV (01 « 1, the second term on the left-hand side may be neglected, 
and we have approximately 

r;'(0 » f J AC sin qo(0 - A sin [ Vo (t) - 0(0]) dt, (20) 

Jo 

where we assume the initial conditions are accounted for in rj o (0- Vo(t) 
varies slowly compared to 0(0, so we treat it as a constant in the integral. 
Our correction term is 

17,(0 = A sin r/o(0 

• f [C - cos 0(0] dt + A cos T7o(0 f sin 0(0 dt. (21) 

Jo Jo 

Equation (21) explains the behavior of the "wiggles" in Figs. 1 and 2. 
If rectangular pulses are used, cos 0(0 = C, and 

11,(0 = A cos j?o(0 / sin 0(0 dt. (22) 

Jo 

For a square wave of period 4T and amplitude 3tt/8, 77,(0 is a saw-tooth 
wave with amplitude? AT cos rj,,(0 sin 3tt/8 = 0.924 AT cos 170 (0- The 
amplitude of the saw-tooth wave is seen to vary as the cosine of 770 (0, 
and reach a maximum amplitude of 0.924 AT, which is 0.023 X tt/2 for 
Fig. 1. This behavior agrees with our earlier observations. The accuracy 
of our approximation in this case is quite good. A plot of r/ o (0 + *?i(0 
superimposed on the figure can not be distinguished by eye from the 
true curve 77(0. 
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We remark that the first term of eq. (21), which is zero for a rec- 
tangular pulse shape, may in general be ignored relative to the second 
term. The first term equals zero at the beginning and end of each pulse, 
and never achieves magnitude greater than AT sin tj (0- The second 
term, however, increases monotonically during a positive pulse and 
decreases monotonically during a negative pulse. The magnitude of this 
term can reach a maximum of § NATa cos r) (t), where 

* = \j Binp(t)dt. (23) 

(The factor f arises because of the possibility of having a message block 
consisting of N +l's followed by a block consisting of N/2 +1'b and 
N/2 — l's). For the case of interest, y (0 will be near zero and N will be 
greater than about 10 or 20. This means that the second term of eq. (21) 
will predominate, and thus eq. (22) may be used for arbitrary pulse 
shapes. (Table I lists the value of a far three important pulse shapes.) 

TV. THE EFFECT OF NOISE 

We now consider a system where both noise and modulation are 
present. We saw in eq. (17) that the zero-order effect of the modulation 
was simply to reduce the locking bandwidth by a factor C. 

Accordingly, we take as our zero-order approximation the solution to 
the stochastic differential equation 

^ + AC sin „ (fl = co - co + j n'{t). (24) 

In this case, we are interested in the steady-state probability density for 
Vo , p(vo)- By comparison with eqs. (6) and (8), the solution can be 
written down immediately. 

where a] , the effective signal-to-noise ratio in the presence of modulation, 
is higher than the zero-modulation signal-to-noise ratio by a factor of 
1/C. 

«• - f - nVTc » 

The average phase shift due to frequency offset is increased. 

. _! co — to 0)0 — co .„-<. 

^ = 8m -Kc~~-Ic- = c- (27) 
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To get a correction term in our solution, we proceed as before, by 
taking the difference between eqs. (24) and (5). We notice that the 
noise term cancels, and we again obtain eq. (18) and the approximate 
solution, eq. (22). For the case of large signal-to-noise ratio, a e , tj will, 
with high probability, be in the vicinity of zero (assuming /3 e ~ 0). Since 
Tj effects 7ji only as the cosine, rj, is essentially independent of t/ in this 
case and is equal to 

tn(t) = A f sin d(t) dt. (28) 

Jo 

rjtit) depends on the particular digital signal being transmitted. 
However, the use of the encoder described earlier insures that 

|*(<)| < § NAT*. (29) 

Thus if each of the components is small, the output phase is seen to 
consist of the sum of three essentially independent parts: 

(i) A constant (o» — o>)/AC resulting from the difference between 

the carrier frequency and the natural oscillator frequency. 
(it) A time varying part which depends upon the digital modulation, 

and which has a maximum magnitude of § NAT<x. 
(Hi) A random part, due to noise, t he dist ribution of which has a 
standard deviation of \Ja, — \/N AC/2A. 

Thus we can write the probability distribution of the reference phase 
approximately as 

pw = e -i5Sp (30) 

where \p, the nonrandom portion of the phase, has a magnitude less than 
or equal to f NATa + |(w - u)/AC\. 

V. CALCULATION OF ERROR PROBABILITY 

If 7) has a known value, the probability of a decoding error, assuming 
equal likelihood detection, is 

P t = \ erfc (p cos 77), (31) 

where p 2 is the signal-to-noise ratio in the bandwidth of the signal. 10 If 
the receiving filter has a bandwidth 2W, then 

A 2 
p2 = NOW' (32) 
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where A is the rms signal amplitude, and NJ2 is the double-sided 
spectral density of the noise. The average error probability P. , is 
obtained by averaging the quantity in eq. (31) over the possible values 
of 77. 

P e = j_ r \ erfc (p cos ,) e ^^ d v . (33) 

This integral was performed using an expansion technique similar to 
that described in Ref. 10. 

VI. DEMONSTRATION OF RESULTS 

In order to reduce the noise as much as possible, the bandwidth of 
the receiving filter in a PSK system is usually set at the value which 
allows the signal to pass essentially undistorted. This bandwidth 2W 
is roughly given by 

m «^?- (34) 

Thus the signal-to-noise ratio in the bandwidth of the loop, a* = 
(4A 2 /N a AC), is related to the signal-to-noise ratio of the received 
signal, p 2 = (A 2 /N 2W), as 

6.4 



«. = 



atc 



(35) 



The effect of noise on the output phase is thus reduced by a factor of 
ATC/6.4 from its input value. We have already seen that the effect of 
the phase modulation on the output phase was proportional to ATo-. 
Thus the size of the quantity AT is important in determining the per- 
formance of the system: it should be kept as small as possible. The extent 
to which this can be done, however, is limited by the need to keep the 
difference between the carrier frequency and the natural oscillator 
frequency small relative to AC. This frequency difference may be re- 
duced by the use of a negative feedback loop. 11 A reasonable value of AT 
presently obtainable in the laboratory for which these conditions can be 
satisfied is AT « 10" 3 . 

Using the value AT « 10" 3 , and the relationship of eq. (35), the 
methods of the preceding section were employed to compute the average 
error probability. Under the assumption that the nonrandom portion of 
the output phase is 10 degrees, the results of this computation for a 
raised cosine pulse shape are plotted in Fig. 3. Also plotted in the figure 
are the curves representing true coherent detection, P e = \ erfc (p), and 
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Fig. 3 — Probability of error as a function of signal-to-noise ratio. 
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differential detection P e = fe~ p \" As can be seen, the curve comes 
quite close to the ideal of coherent detection. The curve corresponding 
to a nonrandom phase shift of 5 degrees was also computed but is not 
plotted in Fig. 3 because it comes so close to the coherent curve that the 
two cannot be distinguished on the scale of the drawing. 

By taking the nonrandom part of the output phase to be zero, it is 
possible to determine the increase in error probability (over the coherent 
case) which is due to noise. We found the error probability to be virtually 
identical to coherent detection in this case. This indicates that the 
effect of the noise on the output phase shift is negligible. 

For reasonably small values of AT(^10 -3 ), and for the range of 
signal-to-noise ratios usually of interest ( > 5 db) the output phase shift 
resulting from noise may be safely ignored: the increase in error prob- 
ability of the proposed system over coherent detection is almost entirely 
due to the effects of modulation and the offset in the carrier frequency. 
Consequently, under these conditions an approximate expression for the 
error probability, which is very nearly correct is P e = § erfc (p cos ^), 
where ^ is the phase shift resulting from the modulation and carrier 
offset. In Fig. 3, for example, the plotted curve is almost identical to 
\ erfc (p cos 10 degrees). 

For AT = 10~ 3 , a total of 289 consecutive raised-cosine pulses of 
maximum amplitude tt/2 with the same polarity are needed to shift the 
output phase by 10 degrees. For AT = 10~ 2 , this number is reduced to 
29. For positive sine pulses of maximum amplitude tt/2, the correspond- 
ing numbers are 241 and 24 respectively. We remark that increasing AT 
from 10 -3 to 10" 2 also increases the effect of the noise on the output 
phase, but that this effect remains negligible. 

In order to demonstrate the effect of the noise on the output phase, we 
must consider an extremely high-noise example. Figure 4 plots the 
average error probability versus the signal-to-noise ratio over a range of 
from —7 to +7 db, for the case AT = 10~\ Nonrandom phase shifts 
of zero and 10 degrees were assumed respectively. As can be seen, the 
zero-phase shift curve almost coincides with the coherent curve for 
p ^ 4 db, even for this large value of AT. 



VII. CONCLUSIONS 

We have analyzed the proposed system of PSK detection for the case 
of random modulation and additive Gaussian noise. If the modulation 
rate is much greater than the bandwidth of the oscillator (AT « 1), 
and if a suitable encoder and decoder are used, we have shown that the 
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Fig. 4 — Probability of error as a function of signal-to-noise ratio. 



system will perform as almost a perfect coherent detector. For a binary 
system, this means a power savings of about \ db over the presently 
employed method of differential detection. This is not very great. 
However, with only a slightly more complex encoder and decoder, this 
same technique may be utilized for higher level systems. (The analysis 
requires only minor modifications.) For a 4-level system, for example, 
the power savings over differential detection is about 3 db which is 
significant. 
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Construction of an experimental encoder and decoder is presently 
being carried out in the laboratory, and tests of the system are planned 
to confirm the theoretical results. 

VIII. ACKNOWLEDGMENTS 

I would like to thank V. K. Prabhu for considerable help and en- 
couragement during the course of this work. The problem was suggested 
by C. L. Ruthroff, who also suggested the general method of approach 
taken in this paper. 

REFERENCES 

1. Hancock, J. C, and Lucky, It. W., "Performance of Combined Amplitude- 

and Phase-Modulated Communication Systems," IRE Trans. Commun. 
Syst., CS-8 (December 1960), pp. 232-237. 

2. Glance, B., "Digital Phase Demodulator," B.S.T.J., 50, No. 3 (March 1971), 

pp. 933-949. 

3. Bussgang and Leiter, "Phase Shift Keying with a Transmitted Reference," 

IEEE Trans. COMTECH (February 1966), pp. 14-22. 

4. Lindsey, "Phase-Shift-Keyed Signal Detection with Noisy Reference Signals," 

IEEE Trans. AES (July 1966), pp. 393-401. 

5. Osborne, T. L., "Amplitude Behavior of Injection-Locked Oscillators," un- 

published work. 

6. Adler, R., "A Study of Locking Phenomena in Oscillators," Proc. IRE, 34, 

No. 6 (June 1946), pp. 351-357. 

7. Ruthroff, C. L., and Stover, H. L., "The Similarity of the Phase-Locked Loop 

and the Injection-Locked Oscillator," unpublished work. 

8. Viterbi, A. J., Principles of Coherent Communication, New York: McGraw-Hill 

Book Company, 1966. 

9. Ruthroff, C. L., private communication. 

10. Prabhu, V. K., "Error Rate Considerations for Coherent Phase-Shift Keyed 

Systems with Co-Channel Interference," B.S.T.J., 48, No. 3 (March 1969), 
pp. 743-767. 

11. Ruthroff, C. L., "Injection-Locked Oscillator FM Receiver Analysis," B.S.T.J., 

47, No. 8 (October 1968), pp. 1653-1661. 

12. Cahn, C. R., "Performance of Digital Phase-Modulation Communication 

Systems," IRE Trans, on Commun. Syst., CS-7, No. 1 (May 1959), pp. 3-6. 



